function exact_approximate_picture(X1, YMatrix1, YMatrix2, YMatrix3, YMatrix4)
%CREATEFIGURE4(X1, YMatrix1, YMatrix2, YMatrix3, YMatrix4, YMatrix5, YMatrix6, YMatrix7, YMatrix8)
%  X1:  vector of x data
%  YMATRIX1:  matrix of y data
%  YMATRIX2:  matrix of y data
%  YMATRIX3:  matrix of y data
%  YMATRIX4:  matrix of y data
%  YMATRIX5:  matrix of y data
%  YMATRIX6:  matrix of y data
%  YMATRIX7:  matrix of y data
%  YMATRIX8:  matrix of y data

%  Auto-generated by MATLAB on 19-Sep-2018 11:36:10

% Create figure
figure1 = figure;

% Create subplot
subplot1 = subplot(2,2,1,'Parent',figure1);
hold(subplot1,'on');

% Create multiple lines using matrix input to plot
plot1 = plot(X1(1:end),YMatrix1,'Parent',subplot1,'Color',[0 0 0]);
set(plot1(1),'DisplayName','var_{dc}/var_{dy}');
set(plot1(2),'DisplayName','\beta_{dcdy}','LineStyle','--');
set(plot1(3),'DisplayName','var_{dc|dy>0} - var_{cd | dy \leq 0}','LineStyle','-.');
set(plot1(4),'DisplayName','\beta_{dcdy>0} - \beta_{dcdy \leq 0}','LineStyle',':');
% Create ylabel
ylabel({'Difference in four key moments';'Exact-approximate solution'});
% Create xlabel
xlabel({'\delta',''});
% Create title
title('CD model: n=3');
ylim(subplot1,[-0.1 0.25]);

box(subplot1,'on');
% legend1 = legend(subplot1,'show');
% set(legend1,'Orientation','vertical');

% Create subplot
subplot2 = subplot(2,2,2,'Parent',figure1);
hold(subplot2,'on');

plot2 = plot(X1,YMatrix2,'Parent',subplot2,'Color',[0 0 0]);
set(plot2(1),'DisplayName','var_{dc}/var_{dy}');
set(plot2(2),'DisplayName','\beta_{dcdy}','LineStyle','--');
set(plot2(3),'DisplayName','var_{dc|dy>0} - var_{cd | dy \leq 0}','LineStyle','-.');
set(plot2(4),'DisplayName','\beta_{dcdy>0} - \beta_{dcdy \leq 0}','LineStyle',':');
% Create ylabel
ylabel({'Difference in four key moments';'Exact-approximate solution'});
ylim(subplot2,[-0.1 0.25]);

% Create xlabel
xlabel({'\delta',''});
% Create title
title('ID model: n=3');
box(subplot2,'on');

% legend2 = legend(subplot2,'show');
% set(legend2,'Orientation','vertical');

% Create subplot
subplot3 = subplot(2,2,3,'Parent',figure1);
hold(subplot3,'on');
plot3 = plot(X1(3:end),YMatrix3,'Parent',subplot3,'Color',[0 0 0]);
set(plot3(1),'DisplayName','var_{dc}/var_{dy}');
set(plot3(2),'DisplayName','\beta_{dcdy}','LineStyle','--');
set(plot3(3),'DisplayName','var_{dc|dy>0} - var_{cd | dy \leq 0}','LineStyle','-.');
set(plot3(4),'DisplayName','\beta_{dcdy>0} - \beta_{dcdy \leq 0}','LineStyle',':');
% Create ylabel
ylabel({'Difference in four key moments';'Exact-approximate solution'});
ylim(subplot3,[-0.1 0.25]);

% Create xlabel
xlabel({'\delta',''});
% Create title
title('CD model: n=4');
box(subplot3,'on');

% legend3 = legend(subplot3,'show');
% set(legend3,'Orientation','vertical');

% Create subplot
subplot4 = subplot(2,2,4,'Parent',figure1);
hold(subplot4,'on');
plot4 = plot(X1,YMatrix4,'Parent',subplot4,'Color',[0 0 0]);
set(plot4(1),'DisplayName','var_{dc}/var_{dy}');
set(plot4(2),'DisplayName','\beta_{dcdy}','LineStyle','--');
set(plot4(3),'DisplayName','var_{dc|dy>0} - var_{cd | dy \leq 0}','LineStyle','-.');
set(plot4(4),'DisplayName','\beta_{dcdy>0} - \beta_{dcdy \leq 0}','LineStyle',':');
% Create ylabel
ylabel({'Difference in four key moments';'Exact-approximate solution'});
% Create xlabel
ylim(subplot4,[-0.1 0.25]);

xlabel({'\delta',''});
% Create title
title('ID model: n=4');
box(subplot4,'on');

legend1 = legend(subplot4,'show');
set(legend1,...
    'Position',[0.755208216708466 0.279367996939038 0.139062497043051 0.143625187983711]);

